UND-HEP-11-BIG04 



On CP Asymmetries in Two-, Three- and Four-Body D Decays 



Ikaros I. Bigi and Ayan Paul 

Department of Physics, University of Notre Dame du Lac 
Notre Dame, IN 46556, USA 



Abstract 

Indirect and direct CP violations have been established in Kl and Bd decays. They have been found in two- 
body decay channels - with the exception of Kl 7r+7r~e+e~ transitions. Evidence for direct CP asymmetry has 
just appeared in LHCb data on Acp(I?° —J' K^K~) — Acp{D^ — >■ tt+tt") with 3.5 cr significance. Manifestations 
of New Dynamics (ND) can appear in CP asymmetries just below experimental bounds. Wc discuss D^^y / 
and Dl/ Ds transitions to 2-, 3- and 4-body final states with a comment on predictions for inclusive vs. exclusive 
CP asymmetries. In particular we discuss T asymmetries in Z) — >■ hih2l'^l~ in analogy with — >■ TT~^iT~e~^e~ 
transitions due to interference between Ml. internal bremsstrahlung and possible El amplitudes. Such an effect 
depends on the strength of CP violation originating from the ND - as discussed here for Little Higgs Models with T 
parity and non-minimal Higgs sectors - but also in the interferences between these amplitudes even in the Standard 
Model. More general lessons can be learnt for T asymmetries in non-leptonic D decays like D — hih2hj,hi. Such 
manifestations of ND can be tested at LHCb and other Super-Flavour Factories like the projects at KEK near 
Tokyo and at Tor Vergata/Frascati near Rome. 



Dedicated to a great physicist and a wonderful friend - Bruce Winstein (1943-2011) 
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1 Introduction 



1 Introduction 

The Cabibbo-Kobayashi-Maskawa theory (CKM) has been shown to represent at least the dominant source of CP 
violation in Bd and neutral K transitions - and that it is insignificant for the matter-antimatter asymmetry in our 
Universe. That asymmetry tells us New Dynamics (ND) with CP violation has to exist. The question is where can 
we find manifestation of such ND and at which level. 

The usual candidate for finding the dynamics underlying the Universe's baryonic asymmetry is neutrino oscil- 
lation with CP violation - and we consider it a very good argument. However we find it wrong for our 'roulette 
game' to place our 'chips' on just one field. Our philosophy for probing CP invariancc for ND is to cover all 
heavy flavour sectors - strange, charm, beauty, top etc. - and also search for EDMs. For charm decays and top 
productions and decays one has hardly any CP asymmetries from the Standard Model (SM). Dynamics of these 
quarks can therefore be well probed for the existence and analyses of ND without too much SM 'background'. In 
charm transitions (unlike for top quarks) one has to deal with the impact of non-pcrturbative strong forces; while 
their impact creates non-trivial 'homework', we want to show that the complexities caused by them offer awards 
too. 

The phenomenology of charm hadrons is richer than for strange ones. There are many more three- and four-body 
final states, and they cover much larger phase spaces. That means that branching ratios are smaller, expansions 
from cliiral dynamics are larger and therefore less reliable than for kaons and hyperons. 

Furthermore long distance dynamics provide larger contributions, over which we have much less control, to decay 
rates than ND. Yet long distance dynamics cannot generate CP violation by themselves. Therefore we conclude 
that rare charm decays can hardly show clear manifestations of ND - unless one observes CP violations, since the 
SM can produce only small asymmetries. 

The rich phenomenology of decays of charm hadrons represent great offers, not just vices: the much larger 
numbers of channels allow more tests of our description of the impact of hadronization. The treating of three- and 
four-body transitions requires much more 'start-up' work; however, once this homework has been done satisfactorily, 
it not only helps in understanding whether ND is behind the observed effects, but also about their features. 

For (and Ac etc.) one can get only direct CP violation and, within the SM, only for singly Cabibbo 

suppressed non-lcptonic decays - except for final states with Ks and K^: for CP odd [even] components in Ks 
\Kl\ generate known CP asymmetry due to indirect CP violation in — oscillation. ND could in principle 
generate direct CP violation in singly- and doubly-Cabibbo suppressed (DCS) channels. The observed Z?" — Z)° 
oscillations can generate indirect CP violation, and many classes of ND can do it well above the SM level. 

In the following narrative, we will discuss CP violation in general for D — > /ii/i2, hih2h^ in Sect. 2 followed by 
the time evolution of the neutral D eigenstates in Sect. 3. Comments on the implication of assuming no direct CP 
violation on measured parameters are made in Sect. 4 followed by a discussion on CP violation in Z?^^ in Sect. 5. We 
look at D'^ — )• h^h^l^l^ in Sect. 6. Using B decays as charm factories is commented on in Sect. 7. The correlation 
of final states in an e+e^ machine is discussed in Sect. 8. then we describe Dl — h^h~l^l~ transitions in relatively 
close analogy to K ^ 7r+7r^e+e^ in Sect. 9 and we comment on other four body final states with a lepton pairs 
in Sect. 10. We discuss D — >■ hlv, hvv in Sect. 11 and comment on D — > h\h2h'ih^ in Sect. 12; then we consider the 
predictions from the class of Little Higgs Models with T parity (LHT) in Sect. 13 in detail and, in general, ND from 
non-minimal Higgs sectors in Sect. 14 before concluding with a summary of our results in Sect. 15. 

2 General Comments on CP Violations in D — >^ '^1^2? ^\^ih-i 
2.1 CPT Constraints 

CPT symmetry tells us the equalities of several classes of transitions; i.e., for singly- and doubly-Cabibbo suppressed 
one: 



V{cJu. qiqiq2q2) = r(cM 91919292) (1) 

T{cu — > dsqq) = T{cu — >■ ^599) (2) 

T{cd — > udqq) = T{cd — >■ udqq) (3) 

T{cd — >■ usqq) ==■ T{cd — >■ usqq) (4) 
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2 General Comments on CP Violations in D — )■ h^h^, hxh^h^ 



r(cs — ?> usqq) = T(cs — > usqq) (5) 

r(cs — usds) = r(cs — usds) (6) 

with q = u,d, s. However, one has to measure transitions of hadrons: 

r{D+/D° -^[S = 0]) = T{D-/D° -^[S = 0]) (7) 

T{D+/D" ^[S = 1]) = T{D-/D" ^[S = -1]) (8) 

TiDt^[S=l]) = r{D-^[S^-l]) (9) 

TiDt^[S = 2]) = r{D-^[S = -2]) (10) 



Final states with S ~ are given with tt's, yy^'^'Tr's, KKn^s, 5* = 1 with A'+tt's, K~^r], . . . and S = 2 with A'+if +7r's 
and K^rjTT. 

Furthermore using strong FSI (ignoring SU(2) breaking in QCD and QED) these relations hold for final states 
separately depending on their isospin /. The equalities of inclusive rates allow asymmetries in exclusive rates, even 
sizable differences. These constraints show some practical usages in 2h, 3h, 4h decays, as we discuss below 

using also G parity. 

For direct CP violation in rate T{D /) vs. r{D /) one needs two amplitudes with differences in both weak 
and strong phases. 

• If the final states are produced by a single isospin amplitude, there can be no CP asymmetry - no matter 
what weak dynamics can contribute whether it be from the SM or from ND. 

• If two different isospin amplitudes can contribute, CP asymmetry can occur depending on SM and the impact 
of ND. In general one expects different landscapes for singly- and doubly-Cabibbo suppressed decays. For the 
leading source of indirect CP violation in beauty transitions one can use the Wolfenstein parametrization. 
However for CP asymmetries in charm decays one has to go to the CKM matrix through 0{X^) and understand 
the differences between other parametrization as emphasized in Ref. [1], where SM does not generate any CP 
asymmetries for doubly-Cabibbo suppressed decays and in singly-suppressed generates ones of order A^. 

2.2 Lessons from LHCb Data on K+K-fn+n- 

In the decays I?" K^K~ vs. K^K^ and tt+tt^ vs. Z)° tt+tt^, one can observe both direct and 

indirect CP violation. The very recent analysis presented by LHCb has searched for CP asymmetry in (mostly) 
time integrated rates in the difference of — >• K~ and 2?° — >■ tt+tt"; its source is (mostly) direct CP violation. 
They find [2]: 

AAcp = Acp{D^ K+K^) - Acp{D^ tt+tt") = -0.82 ± 0.21(stat.) ± 0.11(syst.) (11) 

with 3.5ct deviation from zero. This is the first significant evidence for CP violation in AC ^ dynamics, and it 
is important to infer whether it is due to SM or ND - but has not been established experimentally yet. The next 
topic is more challenging, and just time will not solve it: Can we decide whether we observe the impact of ND? 
There one needs help from theorists - and it will take time-consuming efforts in several ways. 

• To get a realistic SM prediction for CP asymmetries in 2h one has to understand the impact of final 
state interactions. Just diagrams showing them does not mean we can evaluate their impact. Those do not 
tell readers whether they are given by local or even short-distance dynamics or could reflect long-distance 
dynamics. 

Controlling the impact of flnal state interactions is not the strong point of lattice QCD studies. On the other 
hand two hadrons like pions or kaons are beyond the range of chiral dynamics. 

• Looking at diagrams one can guestimate the scale for inclusive direct CP violation within SM of the order 
of 10~^. FSI can change CP asymmetries for exclusive channels significantly. Then one would want to probe 
three-, four-(and more) bodies for larger and smaller asymmetries. 

• If direct CP violation in D'^ — > K^K^, tt+tt^ has been established - and in particular on the level of 
severalxlO"'^ - one thinks about local operators that can be the gate to ND. Such ND would affect CP 
asymmetries in £>+ — AT+ATg (not in £>+ — >■ tt+tt^ as explained below) and in Df — >■ K^(j), K^tt^, Kstt'^- 



2 



3 Time evolution of the neutral D cigenstates 



• It is important to study decays with three-body final states like D — > KKtt, 3tt and Df K~ , 
JC+TT+TT". Analyses of Dalitz plots take not only a lot of experimental work, but also theoretical efforts. The 
'Miranda' procedure [3] can be done in a model-independent way - however one should not ignore theoretical 
help like from analyses of low-energy hh scattering that also include experimental inputs through dispersion 
relations. The calculation effort for a single transition is sizable - and can get help from symmetries arguments. 

• One should analyze DCS decays like Z?" K+tt- , K+tt-ti°, D+ K+tt°, K+tt+tt- and D+ K+Ks, 
Df — s- A"+isr+7r~. The rates for them are obviously lower, the SM 'background' for direct CP violation is 
practically zero, since there is only a single SM amplitude. 

For neutral charm mesons within CKM theory direct CP asymmetry has been predicted around 0(10^"') for 
D'^ — > K^K^ /tt^tt^ with opposite signs (unless FSI affect them significantly) and indirect CP violation, due 
to \q/p\ 7^ 1 and non-zero ^{^Af igD Af), around 0(10"'^); the manifestation of the latter is suppressed by the 
CP insensitive observables xd and yu, which are measured to be around 0.005 - 0.01. The experimental upper 
bounds on indirect CP violation in h'^h^ are well above what the SM can produce. LHT, as a non-ad- 

hoc class of ND, can produce values for indirect CP violation right up to those bounds, while hardly enhancing 
direct CP violation [4]. For direct CP asymmetries in — >■ /tt^h" predictions have mainly been based on 

factorization for (quasi-)two-body final states; looking on diagrams gives only guestimates for the impact of FSI. At 
present one cannot prove that LHCb data are beyond SM dynamics. One has to decrease experimental uncertainty 
- and the theoretical ones too; we will come back to it by emphasizing the importance of the analyses of three- and 
four-body final states. 

The sources of the asymmetry between — > h'^h^ vs. — >■ h^h~ can be differentiated by their dependence 
on the time of decay. To observe indirect CP violation one needs sensitivity to time of decays longer than t^o - 
as illustrated by the CDF analysis [5]. As ND like the LHT models can produce indirect CP violation right up to 
the experimental upper bounds; CDF data has reduced the range of their parameter space sizably. Alternatively, 
models with non-minimal Higgs sectors can enhance both direct and indirect CP violation in these channels which in 
turn makes experimental sensitivity to decay times much larger than tj^o even more important. The flavour-tagging 
of the initial neutral D mesons can be done by the transition of _D*+ — > D'~'tt~^ vs. D*~ — )• Z)°7r~. This allows the 
measurement of the time of flight of the decaying neutral D meson from its birth in general. 

The transitions Z?" — > Kgcf) and l)" — > i^s0 represent a promising place to exhibit a manifestation of ND, since 
the SM cannot generate direct CP violation (apart from Kg ^ K^) [4]. Yet a full analysis has to be performed 
of the Dalitz plot for D°/D° KsK+K'; likewise for ^ Kstitt with inclusion of D° / D° i^s/o(980) 

and D^/D^ — > KsC7{600) contributions; while the total numbers of events in D and D decays depend on their 
productions, the relative ratios in corresponding regions in both Dalitz plots do not. Moreover, direct vs. indirect 
CP asymmetries can be differentiated not only by their time evolution, but also their pattern in the locations in 
the Dalitz plots for three-body final states. 

It is very unlikely that CP violation could produce r{B^ D^D^) ^ T{B^ D'^D^); therefore we can 
assume r{B- -> D°D~) = T{B+ Z)0D+)^ observation of r{B- h+h-D-) ^ T{B+ h+h-D+) can be 
generated by CP violation in h^h^ vs. h^h^ . The same holds true for B^ D^-k^ vs. _B+ D^-k^ . 



3 Time evolution of the neutral D eigenstates 

The decay rates of 1?° and to their respective conjugate states are proportional to: 



\T{D^{t)^,ff = 



\T{D\t)^f)\ 



\AfY 



\Af? ) cosh (y^ — 



\A, 



+25i (J A, ® A*^ sinh {y^±-^ 23 {^Uf ® A*^ sin {x^±-^^ 



1 



-rt 



\A, 



Af\ cosh 



Too 



+2'^{-Af®A*A sinh [yr,— 
V9 V V '^D" 




23 -Af®A*, sin U 



\A, 



Too _ 



td« 



(12) 



(13) 
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3 Time evolution of the neutral D cigenstates 



It is assumed that / and / are final states accessible to both and D'^ . Even the case of Z?" — > — > l'^X~ 
is included, when Ado_^i-x+ = = Aqo^i+x- holds. We use the symbol Af ® A*^ rather than just AjAy. for 
i/iree-body decays the total areas of Dalitz plots depend on the numbers of the produced charm mesons, but not 
the relative ratios between regions in the plot. 

In — > 7r+7r~, K'^K~ , K^tt~ vs. D'~' tt+tt", K'^K~, K'^tt^ the differences between direct and indirect 
CP violations can be shown only due to their time evolutions of the partial widths. For non-leptonic D'^ — > /11/12/13 
vs. -D" — > /ii/i2^3 it was mentioned that the pattern of CP asymmetries can be differentiated by their locations in 
the Dalitz plots. 

Allowing for both direct and indirect CP violation, presently available data lead to [6] ^ : 



D 



= (0.63l°i^) % 



Ar 

yz5 = ^ = (0.75±0.12)% 



2r 



= 0. 



3+0.17 
^-0.15 ' 



D 



= -10.1 



(14) 



i.e., a very different pattern from neutral kaons: 



• The 'mass' eigenstates are close to flavour ones and far from being CP eigenstates, and cuts on lifetimes 
hardly produce the latter. 

• While the experimental value of \q/p\ for D° meson is within one half sigma of the benchmark unity for no 
CP violation, it would also be within one sigma of a value of 0.7 - i.e., large indirect CP asymmetry. 

• ND can produce such large indirect CP violation - and can generate sizable direct CP asymmetries in once 
and twice Cabibbo suppressed transitions, although not necessarily both can be produced by the same class 
of ND. 

• Therefore, one needs excellent time resolution to differentiate indirect and direct CP violation, and one needs 
to compare time evolutions of and transitions. Realistically one can retain only terms through first 
orders in Xj^ and y^,. 



For D° rX vs. D" 



Z+X we can use: 

\T{D^{t) ^ rX)p (X e-" 
\T{D"{t) l+X)\^ cx e-i"* 



1^, 



>i+x\ 



(15) 
(16) 



with |j4j5o^/-_)f I ~ \Ado^i+x\- The indirect CP asymmetry between the 'wrong-charge' semileptonic 
rates are time independent: 



r(ipo -> l-X) - r{D" ^ l+X) _ \q/p\^ - |p/g|2 



T{DO ^ l-X) + T{D" l+X) \q/p\^ + \p/q\^ 
To find the time integrated rates one has to include second orders in Xjj and to get 

^2 



r(D° ^ rx) 



r(£i° ^ l+X) cx 



(17) 

(18) 
(19) 



While those rates are tiny due to the size of x|, 

\qlp\ ^ 1- 



yjj - yet the CP asymmetry might be sizable due to 



^Up to date results can be found in the HFAG website 
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4 Implications of assuming no direct CP violation in charm dynamics. 



Z)" Ks(j) come mostly from Cabibbo favoured transition, where contributions from ND and direct 
CP violation is very unlikely to occur. Therefore we use use |^DO_j.;fg0| ~ \Ad°^Ks(I>\ ^^'^ the simple 
expressions: 



|T(i?0(i)^X5</')poce-r* 

For indirect CP violation one thus gets 

|rp°(t) ^ A-s0)p - |T(J"(t) -> Ks4>)\' ^ 
\T{DO{t) ^ KsW + \T{D"{t) ^ Kscl^W 



1 



t 

t 

TD 



(Vd cos^^ -x^ sin0^) 
(Vd cos <j>jj ~xjj sin0£,) 



(20) 
(21) 



1 t 


Vd 


Q 




P 












) 




P 




q 





(22) 

^ 7^ 0; the term proportional 



i.e., another sign of indirect CP violation that also includes 0^ ^ 0: 

* The term proportional to yj-, is driven by \q\ ^ \p\, yet decreased by 
to X]j is driven by 0^ ^ and enhanced by \q\ ^ \p\. 

Once Cabibbo suppressed decays like D — >■ K~^K~ and D — > tt+tt^ gets both indirect and direct CP 
violation already in the CKM theory, however these are very small asymmetries: differs from 

unity by 10"'* or less, |1 — \q/p\\ and 0^ are generated to about 10^^. A class of non-ad-hoc ND like 
LHT can generate |1 — \q/p\ \ up to even 0.3 and |0^| up to 30° - i.e. up to present experimental bounds. 
However, it is possible, but harder to obtain direct CP violation in these decays much larger than what 
comes from CKM theory. On the other hand, CP violating phases in models with non-minimal Higgs 
sectors can enhance direct CP violation in these channel. 

A'^7r+ time dependent CP asymmetry should be 



- For DCS decays like 



much larger due to the interference with AC = 1 dynamics through Ajjo^x+tt 
^D0->/<'-7r+M_D"->A:-7r+- Finally in the SM one gets 



-/Ado. 



*K+7 



and 



1^ 



D»^K+-K- 



l^50->/^-7r+P and \Ad«.^k--k+? — \AQa^K+^- 



as explained in the beginning. For c — > dsu transitions one gets a combination of A/ = & 1 amplitudes - 
i.e., two different strong phases; ND like the exchange of charged Higgs states can generate a weak phase 
different from the CKM theory. Therefore one can get also a direct CP asymmetry due to |y4£)0_^x+7r- P 7^ 

2 



>X-7r+ I 



Resume for _D" vs. transitions 



1. There are three classes of observables sensitive to CP violation: (i) \q/p\ ^ 1, which is purely indirect CP 
violation, (ii) \Af \ ^ \Aj\ is due to direct CP violation and depends on the final states, (iii) The phase 
between q/p and Af /Af due to the interference between AC = 2 and AC = 1 dynamics; direct CP violation 
will in general induce a different phase for two different final states. 



One can study CP asymmetries in _D*+ 



D"tt+ [1-X]tt+ vs. D*- D"tt- [1+X]tt- and D*+ 
D"7r+ [Ks(l)/K+K- /tt+tt- /K+Tr-]TT+ vs. D*- D"tt- [Ks(t)/ K+ R- /ti+it- / R- tt+]tt- . Such com- 
parisons give us more validations for the measurements and more information on the direct and indirect CP 
violations and the features of the underlying ND. As mentioned above one can analyze B~ — D^D'^ vs. 
B+ D^D+. 

3. Searches for time dependent CP asymmetries have to face the challenge of dealing with yjj and x^, of 1 % or 
less. Furthermore for small effects around 0{1Q^^ — IQ^^) one has to worry about fiavour-tagging. 



4 Implications of assuming no direct CP violation in charm dynamics. 

In a previous work [4] we had suggested that the tt+tt^/A'+A'^ data [5] should be interpreted as being driven 

by indirect CP violation and not by direct CP violation if the ND involved belongs to the LHT-like framework. 
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Direct CP Asymmetries in Decays 



The Heavy Flavour Averaging Group (HFAG) [6] has recently made an analysis of the oscillation and CP violation 
parameters assuming that there is no direct CP violation in charm dynamics. Only two of the listed parameters 
suffered additional constraints under such an assumption, namely, \q/p\ and The quoted values are: 



1.02 ± 0.04 , 0c = {-^-^^-I fi) (23) 



Throughout this article we have looked at two classes of ND which have two different contributions to CP 
violation. 

• One class belongs to the LHT-like framework and can enhance indirect CP violation in charm by orders of 
magnitude above the SM estimate. However, this class of model cannot enhance direct CP violation and 
hence all sources of direct CP violation has to come mostly from the SM. 

• The other class belongs to models with extra Higgs sectors and can enhance both direct and indirect CP 
violation over and above what the SM can generate. 

Of course, there are always possibilities that ND manifests itself as a hybrid of these two classes of models. The 
numbers quoted in eq.23 applies to model that belong to the LHT-like framework and we shall use these constraints 
when discussing the effects of ND of this class. However, an analysis of how the limits on \q/p\ and (t>D constrain 
the parameter space of LHT-like models will be left for a future work. 

5 Direct CP Asymmetries in D'^^-^ Decays 

No matter if direct CP asymmetry in D'^^ — > /l+A'~/7r+7r~ is established or not, there are strong reasons for the 
existence of ND and the study of its shape - in particular in DCS I?^^ decays. 

5.1 Singly-Cabibbo Suppressed Decays 

Ignoring isospin breaking in QCD and corrections from QED there can be no CP asymmetry in — > tt^tt*^ even 
with ND, since there is only one 1 = 2 amplitude and therefore no strong phase difference. 

The same holds for £>+ K^Ks vs. — )• K~ Ks (apart from Ks not being a pure CP eigenstate): c — )■ ssu 
represents a A/ = 1/2, A/3 = +1/2 transition and therefore cannot mix with _D+ — >■ 7r+7r° in the limit of isospin 
symmetry. 

For Z)+ — > A'sTT^" vs. D~ — > Kstt~ one can (even in addition to Ks 7^ K^) have a difference due to ND 
affecting the DCS transition c — > dsu that interferes with the SM amplitude c — > sdu for — )- Kstt'^. Due to the 
known CP asymmetry in the neutral K system this asymmetry can be predicted to be: 

— ; ^-i ^ '- ~ 2Rc(eK) ^ 3.3 x 10"^ (24) 

r{D+ KsTr+) + T{D- Kstt-) 

If ND contributes significantly to the DCS amplitude c dsu, it will probably change this prediction; sizable FSI 
between Ar°7r+ and A'°7r+ are expected. 

In analogy to decays one predicts a CP asymmetry in Df — > KsK^ vs. — > KgK^ of 3.3 x 10^^; 
an impact of ND on the amplitude for c — >■ dsu can affect this prediction. For Cabibbo suppressed transitions 
£)+ AT+tt" vs. Dj — )- K^tt'^ one predicts asymmetry of order 10~^ to lO^"' in the CKM theory; for 1)+ — > tt+ATs 
vs. T^^Ks one expects asymmetry up to 4.3 x 10^'^ due to direct CP violation in AC = 1 CKM dynamics. 

There are subtle features here: 

• Through order there is no weak phase for these DCS transitions within the SM; therefore there can be no 
CP asymmetry on the quark level, never mind FSI. 

• There are weak phases of order A"' for Cabibbo favoured transitions. However the same weak phase controls 
the amplitudes, therefore no CP asymmetry can be produced at the quark level. 

• Yet one measures transitions at the hadronic level: Z?+ K^h's and Z?+ — > K^K^h's (with h denoting 
it's and 77/77'), but also D+ -> Ks,l + h's and D+ K+Ks,l + h's/Ks,LKs,L + h's. Those mix Cabibbo 
favoured and doubly suppressed amplitudes. 
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6 D° h+h-l+l- vs. D° h+h-l+l- 



5.2 Doubly-Cabibbo Suppressed Decays 

In — > A'^tt" decays the final states contains interference between 7 = 1/2 and / = 3/2 amplitudes with different 
strong phases. In the SM there are no weak phases through order - i.e., no practical CP asymmetry. However 
ND can generate sizable weak phase to compete with doubly-Cabibbo suppressed amplitudes. The landscape for 

is even richer, as discussed below. 
For Dg decays one can measure Df K^Ks - but asymmetries can 'hide' in Cabibbo favoured decays. 
However one can probe £)+ —J- K'^K'^tt^ vs. K~ K~tt^ . 

5.3 D^^^ ^i^2^3 

The decays of charged D mesons to final states with three pscudo-scalars give us much more information about the 
dynamics underlying direct CP violation than studying final states with two pseudo-scalars, where one can measure 
only differences in the partial widths. For three-body final states one can measure - and it has been done - the two- 
dimensional Dalitz plot with the known resonances. There one can differentiate strong vs. weak phases and provide 
more information about the possible impact of and the nature of ND. While performing a full Dalitz plot analysis 
requires more data and time ~ even using the 'Miranda' procedure [3] - it gives us a great deal of information on 
the ND; in particular one has to include the informations about contributions from scalar-pseudoscalar final states 
like Ka, Kk and na. Most CP asymmetries do not depend on the relative vs. 13^^ productions. There are 
three classes of CP asymmetries: 

I. A difference in 13 ^ Rh vs. D Rh; 

II. Difference in the interference between D — > Rihi and D i?2^2 vs. D — > Rihi and D -R2^2; 
III. Difference in true three-body final states or broad resonances. 

For classes I and II one can use Breit-Wigner prescription for resonances even for an accurate analyses; however for 
class HI one needs some input from theory to use dispersion relations from /ii/i2 scattering at low energy. Several 
studies are at work. 

In our previous paper [4] we have shown that LHT models for ND do not have a good chance to manifest 
themselves through direct CP violations in charm decays. However ND models with a charged Higgs sector could 
have a sizable impact, in particular for final states with scalar resonances like Kfo{980), Ka{600), 7rfT(600) etc. 



6 ^ h+h'l+l- vs. h+h-l+l- 

In K — > 7r+7r^7 one has partial widths and the T odd correlation (ty • (p+ x p-)), where describing the photon 
polarization. The most realistic way to infer is based on K'^/K'^ — > tt^tt^j* — s- 7r+7r~e+e^ and in measuring 
the angle $ between the 7r+ — tt~ and e+ — e~ planes. 

The largest CP asymmetry in kaon decays has been found in Kl — > e^e^7r+7r^, namely this T odd correlation, 
a prediction by Sehgal and Wanninger [7] 

^T|thoor.y= (14-3 ±1.3)% (25) 

that was measured later by the KTeV and NA48 experiments 

^tL.p = (13.7 ±1.5)%. (26) 

The large effect is generated by the small quantity \ri^ | in the Kl wave function. As usual a price has to be paid 

for such an enhancement, namely a tiny branching ratio of (3.11 ± 0.19) x 10~^. In the context of beauty physics, 
such kinematically nontrivial CP violation stemming from triple products of spin and momenta was discussed in [8] . 
This can bring about some hope of finding an enhancement in _D — > l'^l~'h^h~ too with I = e, fi and h = K , tt 
again with the normal price of a tiny branching fraction. 

It should be noted that the analogy between such D° with K'^ decays is only qualitative: 
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• The lifetimes of the two mass eigenstates in the — K'^ system are very different so one can separate Kl 
and Ks easily; therefore the notation is obvious. The lifetime difference in the D'^ — complex is only 2% 
or less and hence one has to compare D° and Z)° time evolutions carefully. To identify a T odd correlation as 
a CP asymmetry one has to compare them in and transitions and show that they are not of opposite 
signs. 

• The mass eigenstates Kl and Ks are almost aligned with the odd and even CP eigenstates. For the — D° 
complex both mass eigenstates might contain sizable odd and even components in the presence of ND. 

• Unlike for the neutral kaon system, direct CP asymmetry might not be smaller than the indirect CP effects 
for the neutral D mesons. Like in the past, we shall continue to assume that the contribution from one class 
of ND, namely LHT-like ND, to direct CP violating dynamics is tiny compared to the indirect one. If reality 
proves to be in the contrary, it will only serve to enhance the effects that we speak of while a quite different 
classes of ND may have to be considered when models like the ones with non-minimal Higgs sectors might be 
good candidates. 

• While the indirect CP violation was known from Kl — > tt+tt^ , we have only an upper bound for transitions. 
Furthermore it is (and was) known that CP asymmetries in neutral kaons is dominated by indirect CP 
violation. However we do not know what pattern holds for indirect vs. direct CP violations in neutral charm 
transitions. 

• K decays can be treated in chiral dynamics with 'soft' pions but this can hardly be trusted in D decays to 
pions, kaons and photons. In any case one cannot expects similar effects. 

• At the same time we can see also virtues, not just vices, namely to have D — > Knl'^l" , KKl^l" , ■mrl'^l^ with 
I = e, ^ and D — > ttttj, KKj, Kttj for comparisons of Ml, El and IB amplitudes to check predictions. 

In terms of the flavour eigenstates the decay can be expressed as: 

-^T{D°{t) h+h-l+r) = Tl{t) cos2$ + Vl{t) sin2$ + Tl{t) cos$ sin$ (27) 

-^r(Z)°(t) ^ h+hrl+r) = f ?(<) cos^$ + f ^(i) sin2$ - flit) cos$ sin$ . (28) 

The change in sign of vs. is due to $ — > — $ under a CP transformation as CP symmetry implies Pi 2 = Pi, 2 
and Pa = -fa. 

CP violation, given by Pi. 2 7^ Pi, 2 can be driven by direct CP violating dynamics only. However, ^ Ps, can 
be driven by both direct and indirect CP violating dynamics. The asymmetry can be written as: 



r{D°{t) ^ h+h-l+l-) - T{D\t) ^ h+h-l+l-) ^ ,„,^ 
V{D°{t) ^ h+h-l+l-) + T{D^{t) ^ h+h-l+l-) Too 



(29) 



where a^p is the time independent part which is driven by direct CP violation and a™p is the time dependent part 
which is primarily driven by — oscillations but can have contributions from direct CP violating dynamics 
too. In the absence of sizable direct CP asymmetry in charm dynamics, i.e. \T{D^ — )■ h+h~l+l~)\ = \T{D^ — >■ 
h+h-l+l-)\, it can be safely assumed that the entire asymmetry depicted in Eq.29 is driven by oscillations. However, 
there is some subtlety that needs to be considered here. The time dependent term a™p is given by 

acp - \ [to5R(A - A-i) - Xo^{\ - A-i)] , A = l T{D°^ h+h-l+r) 
CP 2 L-^-^ ^ ^ D \ )\ 1 pT{D^ ^ h+h-l+l-) ^ ' 



In the absence of direct CP violations we have 

A 



g*(0+<5) 



(31) 



The angle cj) is the phase of q/p and comes from mixing of the flavour eigenstates and hence contributes to 
oscillations. However, the phase 5 is final state dependent as it comes from final state interactions (FSI) and is the 
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phase difference between T{D" h+h-l+l-) and T{D" h+h-l+l-) which can be present even in the absence 
of direct CP violation in this decay channcL This argument holds true for any final state if oscillations can drive 
a CP asymmetry in it. Hence a genuine difference in a™p between different decay modes can be an indirect signal 
of direct CP violation. 

While the observation of Z?" — _D" oscillations seems established, the relative size of xd and yo is not clear 
yet. Before these experimental results in 2007, most authors had argued that the SM predicts xd, Ud ^ ^ ^ 10~^ 
— yet not all. In 1998, xd, Ud < 10^^ was listed, admittedly as a conservative SM bound [9], together with 
a question: How can one rule out that the SM can not produce 10~^ 1^ < 10~^ (corresponding to xdi Ud 
~ 10"'^ — 10~^). In 2000 and 2003, an SM prediction obtained from an operator product expansion (OPE) yielded 
Xd, Ud ^ 0(10^'^) [10]; later a more sophisticated OPE analysis was done with similar results [11]. Alternatively 
in 2001 and 2004, an SM prediction on — oscillations was based on ^[/(S) breaking mostly in the phase 
space for yD and then from a dispersion relation for xd [12]. 

Even though the present experimental results on xd and y^ can be accommodated within some available 
theoretical SM estimates, ND can play a significant role here [13-15]. In particular in [13] such has been analyzed 
for LHT-like dynamics, if future data confirm that xd indeed falls in the range of 0.5% and 1%. While SM long 
distance dynamics could accommodate a value in that range, it could not be ruled out that ND could contribute 
half or quarter of the value of Xjj considering the reasonable theoretical uncertainty in the predictions from long 
distance dynamics. Yet indirect CP violation in — oscillations would provide us a clear signature for the 
manifestation of ND. 

From the theoretical SM estimates for xjj and yjj , which arc hardly more than guesses since the SM predictions 
depend strongly on our theoretical treatment of LD dynamics for xd and yu (or its lack thereof) and our treatment 
of the extraction of the CKM phase as it enters charm decays [10, 15] one gets a SM "prediction" of: 

«b"^(i)lsM~O(10-') , (32) 

while the same for direct CP violation is most likely 

acpWlsM^O(10-^) (33) 

Xd = 0.63%, Yd = 0.75% 




0.6 0.7 0.8 0.9 1.0 1.1 1.2 0.6 0.7 0.8 0.9 1.0 1.1 1.2 



Iq/Pl Iq/Pl 

Figure 1: Contour plot of the indirect CP asymmetry parameter a™p . The plot on the left allows for any arbitrary 
value of S in the range {— tt, tt}. The ellipse marks the possible asymmetry assuming no direct CP violation with 
the dot marking the central values under such an assumption. The \q/p\ and (j) values used here are the 95% CL 
limits as quoted in [6]. 
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for the final states where h+h' = -k+tt- /K+K~ . If we consider D° K+n-l+l- vs. D° K'tt+I+I-, SM 
cannot generate direct CP violation in this DCS channel. The contribution to the decay rates from the SM are 
smaller too and hence there can be possible large ND intervention here. 

In Fig.l we show how large a^p can be for the current central value of and yjj. We have varied q/p and 
(j) over the current 95% CL values. For the figure on the left we have allowed an arbitrary variation of S over the 
entire range of {— tt, tt}. In Fig. 2 the same is depicted but over variations of xd and t/d over their current 95% CL 
values with a choice of q/p and (/) where the a™p is maximized for these parameters. It can be clearly seen that at 
most 

a^^ ~ 10-2 (34) 

Such time dependent CP asymmetry can be measured in hadron machines as has been done for — > tt+tt^, K^K~ 
at CDF recently [5] . What is important to distinguish between indirect and direct CP violations is that experimental 
analysis of such time evolution should be conducted over large values of (t) / t^o . Possibilities of measuring the 
asymmetries at varying values of {t)/T]jo ranging from {t)/T£)a = 1 to {t)/T£,o ^ 1 should be looked into. Later we 
will show that ND can come close to such upper bound, namely from LHT-likc models or models with non-minimal 
Higgs sectors. However, assuming no sources of direct CP violation in charm dynamics puts additional constraints 
on models within the LHT-like framework. 

Akin to the analysis done for — >■ 7r+7r~e+e- in [16] a "time-resolved" analysis can be done for the system 
too. The charm specific analysis would display the evolution of the asymmetry from its zero value at short time 
scales to its maximum value at large time scales. A spectrum integrated asymmetry would display the interference 
between the two mass eigenstates and an asymmetry between the decay rates of vs. can also be analyzed. 
At this time such a study is not feasible as we have very little knowledge of the strong phases involved in the decays 

— > h^h^. However, such an analysis would be quite instrumental in determining what scales of (t)/r£)0 to 
probe for experimentally. It should be kept in mind that due to the tiny lifetime difference between the two mass 
eigenstates of the meson, time resolution of the order of the difference might be necessary to do a time resolved 
study, a capability that shall, hopefully, exist in the future. 

The plots clearly show that the observable CP asymmetries cannot be generated above few x 10"'^ in general. 
That is not surprising, since the CP insensitive parameters xd and yo - where ever they come from - are around 
0.005 - 0.01: ND cannot generate more than that range. 
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7 Decays as a Factory of Neutral D Mesons 

The decays of charged B mesons can act as a 'filtered' factory for neutral D mesons with sizable branching ratios: 

BR{B+ D°D+) = 0.010 ± 0.0017 (35) 
BR(B+ ^ D°D*+) = 0.0076 ± 0.0016 (36) 

Searches for small CP asymmetries in neutral D transitions in three- and four-body final states should be possible 
at future Super-B factories, where at least partially time integrated asymmetries in final state distributions should 
be possible. 



8 On e+e- D1D2 

If one would have correlated — pairs, one could find larger effects. It can happen in e+ — annihilation 
realistically at the '0"(377O) resonance. 

Possible CP asymmetries in D transitions are very small and the distinction of indirect vs. direct CP violation 
over time evolutions very subtle. Then the question arises that quantum correlations in e'^e" — > D^D'^ — > DsD^ 
might help us here in some qualitative analogy with e+e~ ^(1020) K^K^: there one can filter the mass 
eigenstates using e+e~ — )■ KsKl [tt'^tt~]ksKl as a flavour-tagged production of Kl — >■ 27r, e+e~7r"'"7r~ due to 
tks < '^Kl ■ For 

e+e- ^ ^"(3770) ^ D" D+D^ / D1D2 ^ fafb (37) 

one finds the DD in a highly correlated C odd P wave configuration. Bose-Einstein statistics forbids the transfor- 
mation D'^D'^ D^D^ /D^D^ although correlated oscillations are possible. Therefore one wants to measure final 
states fa and that can be fed by both coherent sum of & D'^ or & £)_ or Di Sz D2: 



BR(i//'(3770) ^ D''DyD+D^/DiD2 ^ fafb) 



BRiD ^ /a)BR(D ^ fb) X 



(2 - xl + vl) \p{fa) - + {xl + yl) 



p 


2 


q 





1--Pifa)-P{fb)' 
P P 



(38) 



This rate is very small due to small values of BR{D fa) and BR{D fb). However there are several reasons 
to consider them: 

• CP is violated by the existence in time integrated rates of final states that both have CP even or odd quantum 
numbers in e+e~ — >■ [fa){fb), see Eq.(38), which provides a distinct statistical advantage over the measurement 
of CP asymmetries as a smaller number of events suffices. 

• One cannot measure some phases in any other way. 

• The rate of e+e" V'"(3770) D^D^ / D+D_/ D1D2 fafb have small backgrounds. 

• The intermediate '0"(377O) state is a CP even state. In principle even a single event where the final state is 
even from the decay of a £)"Z)° pair establishes CP violation. 

• Such CP violation can occur either due to the existence of direct or indirect CP violating dynamics or both. 

• The existence of the transition where both fa and fb are either CP even or CP odd proves the existence of 
CP violation in the dynamics regardless of whether fa and fb are the same or different final states. 

• Indirect CP violation in neutral kaon decays has been manifested in a much larger asymmetry in Kl — >■ 
7r+7r"e+e" decays through a final state distribution of the 7r+ — tt" plane vs. the e+ — e" plane - at the 
price of the much smaller rates! Something might happen for final states distributions in neutral D decays 
in £) — >■ h+h~l+l~ or — >■ if +if "tt+tt" as discussed in the next section. 
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As a consequence of EPR correlations [17], tagging through final states h+h~ ~ like tt+tt" or K+K~ - one gets 
a D- beam with some purity; likewise tagging Z?_ through final states Ks4' one gets a beam. However, the 
correlation between the two eigenstates last only till one of them decays. After that the other eigenstate that has 
not decayed is prone to oscillations and the time evolution of the state has to be carefully studied. If both fa and 
fb are CP even [odd] it can be clearly inferred that, at the time of the decay, it was a D^[D-] which decayed first 
and the correlated state which was D- [D+] oscillated in to the other CP eigenstate of hence displaying the 

existence of indirect CP violating dynamics. Of course, we are assuming here that the effects of direct CP violation 
is tiny in comparison. 

For the case of the DsDl complex produced at threshold, tagging thorough final decay products is no longer 
trivial as in the case of the kaons since the mass eigenstates are no longer CP eigenstates to any good approximation. 
Experimental results allow for 0.8 < \q/p\ < 1 within about a la range. This, however, does not mean that an 
analysis on the basis on mass eigenstates is impossible. A good measure of the CP impurity in the mass eigenstate 
can possibly allow tagging of the mass eigenstate based on a statistical analysis. 

In e+e~ annihilation close to charm threshold one can use correlations between and Ds transitions like 
l+l+X vs. X and l+hih2X vs. l~hih2X without a time filter. The fact that both Dl and Ds can contain 
sizable even and odd CP components leads to complex time evolutions as: 

sin^ (ii - t2) (39) 
sin2 ^T{t^ - t2) (40) 

It is also quite informative to study the time evolution of the states [l^X]t-^ [Ks(l)/h+h~]t^ vs. [l+X]t-^ [Ks4>/h+h^]t2 
or [l^ X]t-^[K+Tr^]t2 vs. [1+ X]t-^[K^TT+]t2 from a pair of correlated D^D^ mesons, as these show a diflFerent depen- 
dence on Xj-,. The asymmetry is proportional to smxjj{ti — t2)/T£)0 as the neutral D meson pair is produced in a 
C odd wave in a e+e~ machine and hence vanishes when integrated over ti and t2. However there is one way out. 
If one studies 

e+e- D^D^j (41) 

then the correlated D^D^ are produced in the C even state. The time dependent asymmetry is now proportional 
to su\Xj-,{ti +t2)/T£,o and hence does not vanish when integrated over ti and t2- The asymmetry is proportional 
to 

2x^3 (^^p{f) 

which is the same as that for incoherently produced neutral D meson. 

The e+e~ machines at the B factories currently run at the T(4S') resonance generating a large number of neutral 
B mesons. As mentioned before, the decays of these B mesons can be used as a "factory" to generate neutral D 
mesons in their fiavour eigenstate using the charge of the associate Ds produced as a flavour tag. In such an 
analysis the time dependent CP asymmetry can be measured although the neutral D mesons that are produced in 
this manner are not very highly boosted, unlike in hadron, machines and hence this allows access to values of only 

t/TDO ^ 1. 

9 h+h-m- 

Tempted by this enhancement in K]^ decays we studied the analogues transitions D — > h+h^l+l^ with h = t:,K 
and I = e,fi. Indirect CP violation in Kl decays can be described by JiK{KL tt+tt~) ~ 2 x 10~^ or r(ii'i — >■ 
TT+TT^)/r{KL — > 7r+7r~) ~ 6 X 10"'^ - i.e., a very small effect (and given by CKM theory at least as the leading 
source). For indirect CP violation in transitions CDF has placed a 95 % bound of 1.4 • 10~^ x {t)/T]jo [5]. At 
first sight that number for CP violation in D'-' mesons looks to be in the same ballpark for the underlying dynamics, 
although CKM theory cannot reach it. 

9.1 The Decay Width of Dl h+h-l+r 

As in the case of kaons, the decays D^ — > h+h^'j can occur through two distinct channels. Firstly, it can occur 
through an internal bremsstrahlung associated with the CP violating decay D^ — >■ h+h~ . It can also have a direct 



-r(ti+t2) 
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emission component, mostly a Ml transition with a possible El component coherent with the bremsstrahlung 
component, both being CP conserving as shown in Eq.42. In e~^e~ — > DD - in particular, close to charm threshold 
- one finds that the final state has to carry even CP parity with very high accuracy due to dominance of one over 
two 7* production. Therefore the pair of neutral D mesons hadronize as DlDs as for e+e^ — > KsKl- First we 
consider 

^ h+h- ^ h+h--f and Dl ^^i^ h+h--f, (42) 

where D/, — >■ h^h~ can be produced from indirect and direct CP violation. The final state contains both a CP 
odd state created by the Ml transition and a CP even state due to the IB component. However, as long a the 
polarization of the final state photon is not observed, this interference cannot be measured [7]. The only way 
to observe this interference is to take the photon off-shell in which case the information about CP and T odd 
polarization of the virtual photon is retained by the final state leptons and the asymmetry manifests itself as an 
angular correlation between the and h~^h~ planes. The photon polarization can be probed best for off-shell 
photons 

Dl h+h--f* h+h-l+r ; (43) 

in addition to the IB, El and Ml amplitudes one also has to include the charge radius term Dl Dsl^l^ — > 
h'^h~l^l~ which contributes to the decay rate but not to the asymmetry. In what follows we will follow the notation 
of Sehgal and Wanninger in [7]. The amplitude for Dl — > h'^h~l'^l~ is described by 



T{Dl h+hrl+l-) = e\T{Ds ^ h+h- 



\ [eP, {P ■ k)k,] i + ^ \{P.k)p+, {p+ ■ k)P,] 



11 11 



^M(fc_)7^i;(fc+) (44) 



where k ^ k-^- + fc_ is the sum of the final state leptons' momenta, p+ and p- the momenta of the /i+ and h 
respectively and P is the momentum of the incoming D meson. The differential decay rate is given by 



' dy / dxAKl,x,y) / — F(x,y,$) 



T{Ds^7T+7r~) 167r2A^ (l,Ai2,M^) 74^= Ji^.^ Jo 27r 

(45) 

Here x — {p+-\-p-Y /m?jj is the normalized invariant mass of the final state mesons h = w/ K , y = (fc+-|-fc~)/m|) 
is the normalized invariant mass of the lepton pair I = e/fi and $ is the angle between the plane containing the 
meson pair and the plane containing the lepton pair. The normalized pion and lepton mass is given by 

mh mi 

^ = , V = (46) 

niD mo 

The function F{x, y, $) can be split into three distinct parts [7]. 

Fi{x, y, $): The terms in this part is proportional to either sin^ $ or cos^ $ and contributes to the decay rate and 
not the asymmetry, hence they contribute only to Fi and F2 in Eq.49. 

F2{x, y, $): The terms in this part contribute to the asymmetry only and not the decay rate and hence contributes 
to F3 only in Eq.49; as they are proportional to sin 2$. 

F3{x, y, $): The part proportional to cos 2$ which contributes to neither of the above. 

In calculating the decay rate, terms proportional to the lepton mass have been ignored which holds good for 
the electrons in the final state. For muons this can introduce an error of at most 0(5%). Our estimate of the 
branching fractions using parameter ranges that have been explained below, put them to 

BR{D TT+TT-l+r) - 10"^ 

BR{D ^ K+K-l+l-) - 10^1° - 10"^ (47) 
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9.2 The CP Asymmetry Parameter in Dl — >■ h~^h Z+Z 

The interference of the CP violating IB and El amphtude with the CP conserving Ml amphtudes wiU yield 
a circularly polarized photon and leads to a triple correlation between the hadrons' momenta and the photon 
polarization 

Pi = ■ iPh+ ^ Ph-)) , (48) 

which is CP odd. The CP violating effect appears as a correlation between the e'*'e~ and h^h^ planes. The 
differential decay rate with respect to $ can be written as: 

—T(Dl h+h-l+r) = Fi cos^$ + r2 sin2$ + Pg cos$ sin*. (49) 
d$ 

It is easy to see that the term proportional to cos $ sin $ changes sign under CP and T; Ps thus represents CP 
violation. Non zero asymmetry is implied by J J J F2{x,y,^)dxdyd^ ^ with the integral in $ done over the 
difference of the unions of opposing quadrants. Thus, the asymmetry parameter is defined as 



A 



^ C -(ri + r2) 



^^■^ (50) 



A detailed look at the form of Fix, y, $) in [7] will make it clear that the asymmetry comes from the interference 
of the CP conserving Ml amplitude and CP violating El and bremsstrahlung components. One should note that 
the CP/T odd correlation is controlled by Vh+f^-, which enters through .ggp, as wc shall see later and is given by 



Vh+h 



^ T{Dl ^ h+h-) 

T{Ds^h+h-) ^ ' 



which can arise from both indirect and direct CP violation; also one should understand the dynamical situation is 
very different from decays: 

• Unlike the known situation in K decays direct and indirect CP violation in D transitions are of similar 
strength within the SM and are weak. In ND they could also be of similar strength and sizable - they could 
reach the upper experimental bounds listed in Eq.(14). In one class of ND, namely LHT, one finds that 
indirect CP violation could be large, yet direct CP asymmetries weak. On the other hand, models with 
multiple Higgs doublets can contribute sizably to both direct and indirect CP asymmetries. 

• One gets more channels where CP/T odd correlations can exist and even measured, namely Dj^ — > tt^^tt^/^/^, 
K~^K^l^l^ with I = e, fjL. One can compare vs. Vk+k- ' ^ difference would show direct CP violation, 
but 11^+^- — Vk+k- 'io^^ prove the absence of direct CP violation. However, one can compare ?7^+^- 
vs. 77^0^0 as usual. One could also study e+e" D°D° Ds + [iCs /o(980)] and probe Ds R'sn^l'^l', 
Ksp°l+l^ for a CP/T odd correlation. 

• While the correlation is controlled by rjf^+i^-, it is also governed by other factors that depend on the impact 
of strong forces. For — ^ 7r7r7(*-' and — ^ Kse^e~ one can apply chiral dynamics with small spaces, but 
there are more problems and uncertainties for D transitions. 

We know that the differences between Z?" h^h^ and h^h~ arc small. Indirect CP asymmetry are 

suppressed by and yjj of about 0.005 - 0.01. In preparing Dl samples, one can find indirect CP violation in 
Dl — >■ h^h~ without this penalty. One might enhance the strength of CP violation in Dl — > h'^h~l'^l~ from ND 
at the price of much smaller branching ratios - if the strong forces collaborating as they did for Kl — > Tr~^ tt~ e'^ e~ . 



9.3 A Note on the Short Distance contribution to — >■ h'^h 

In addition to the contributions discussed in the previous sections, there can be a short distance contribution 
to Dl — >■ h^h~l~^l~ analogous to similar effects in — >■ 7r+7r~e"'"e~ as discussed in [18]. The short distance 
contribution can be classified into three types. There can be vector and an axial current components stemming 
from a local interactions of the current U7^(l — 75)0 with ^7^/ and ^7^75^ respectively. These amplitudes can possibly 
contribute to both the decay width and the CP asymmetry. Within the standard model this contributions will be 
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highly CKM suppressed and a rough estimate clearly shows that it will be several orders of magnitude smaller than 
any of the contributions discussed above. The third component leads to the CP violating El amplitude that has 
been discussed above. This component is also CKM suppressed within the SM. 

ND like LHT can have large effects on these amplitudes and can also lead to enhancements of several orders of 
magnitudes over the SM in certain parts of the parameter space. This is not unusual as such has been seen even for 
other short distance operators studied previously for other decay channels of the neutral charm meson. However, 
even with such enhancements, long distance effects remain the dominant contribution to this decay channel. 



9.4 Estimation of Parameters 

To complete the analysis we have to estimate the magnitude and phases of the following five parameters from 
experimental measurements. 

(i) The decay constant f^: The decay constant fg is defined by 

T{Ds ^ h+h-) = (i (52) 

IbTrmr, \ jtiJj / 

we will only require the phase information of fg . 

(ii) The bremsstrahlung parameter gg^^: is defined as 

9BK^Vh+h-^y arg(gBR) = $±+'5o (53) 

where 

D{h) ^ {h+h-\Hw\DL) _ , ( D{h) \ ^ D(h) 

"^'^^i^- - {h+h~\Hw\Ds) ~ ^ + ^' ^ y'^^'^-J - ± ^ 

Assuming that the magnitude of fy^-f is driven purely by oscillations and the effect of direct CP violation is tiny 
we get 

(55) 

which is an approximation good to about 0(10%) since in the adopted CKM phase convention the assumption that 
T{D° ^ /) ~ T{D" f) holds good [13]. Taking the 95% CL values from HFAG [6], 



e {0.6, 1.28}, e {-27.5°, 8.4°} =^ |??± | e {0, 0.3}, $± e {-TT, tt} (56) 
The phase do in Eq.53 is the hh scattering phase in the I = J = channel at y/s = mo- 



(iii) The Ml parameter g^^^'- A measurement of the direct emission component oi ^ h^h 7 is necessary 
to estimate g^^ 



-{9miY TI^£r U-— =(3mi)'^mi(™/.) (57) 



T{Ds^h+h-) m^&poV TTiD 

where 

Am? \ ^ 0^ 



1-^ 7^ , /3o=/3U=o, iOma^^rriD^ (58) 

\ mf) — 2m£)W / 2 

With a cutoff of ujmin ~ 20Mey, /mi("1/i) is of 0(10^^) for h ~ tx and 0(10^'') for h = K. Hence, although we 
cannot calculate g^^ at this point for the lack of experimental data, we can estimate it to be of 0(1) - O(10~^). 
The phase of g^i is 6i{sh) + 7r/2, where Si{sfi) the hh scattering phase in the vector channel at a cm. energy of 
y/s^ and the phase shift of 7r/2 has to be introduced to allow for the phase convention used for g^^ [7]. 
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9 Dl ^ h+h-l+l- 



(iv) The El parameter g^^: The El component is distinct from the IB component and is a CP even 
component in the — ^ h^hr^ amphtude. The phase of g^^, however, depends on whether it comes from the 
CP impurity in the wavefuntion or from an intrinsic CP violating dynamics in the CP odd part of D^. In the 
former case the phase of g-^^ relative to that of g^^^ is given by 



arg ( ^ U <i>e « 'f^ - - (59) 
\5mi/ ^ 

In the latter case g-^^ is out of phase with g^^^ and hence does not contribute to the asymmetry. As in [7] we will 
continue to use the first scenario and not the second. Also, as in the case of the kaons, it is logical to assume that 
g-^^ will be tiny compared gy^^ and in the lack of any evidence to the contrary we shall maintain the same. It is 
unlikely that ND contributions to the El amplitude can enhance gEi significantly. Hence we shall take 



TT 



5ei 



5m 1 



0.05 (60) 



(v) The charge radius contribution parameter g^: The amplitude connected to the charge radius does 
not contribute to the asymmetry but contributes to the total decay rate. This parameter depends on the estimate 
of the D meson charge radius, which is estimated to be twice as large as the neutral Kaon charge radius. Along 
with the fact that mo ~ 3.6 x tuk, we estimate to be 

9p = -\{R^)ml^7.S (61) 

The phase of gp is (5o(s/i) which is the hh scattering phase in the / = channel at cm energies of ^/Wh■ It should be 
noted that we do not consider any modifications to gp due to the off-shell behavior of the Ds — >■ h^h~ amplitude 
as pointed out in [18] 

9.5 Numerical Results 

From the previous section it should be quite evident that a significant degree of uncertainty plagues the estimation 
of the required parameters. An exact (or near exact) calculation of this asymmetry requires far more experimental 
input than we currently have. The largest uncertainty, by far, comes from the estimation of the charge radius 
parameter g^ which is only a theoretical estimate of a purely long distance contribution and hence way beyond our 
theoretical control. There are also large uncertainties in the estimation of both the magnitude and the phase of 
5bri the constant parametrizing the bremsstrahlung contribution to the amplitude. The uncertainty stems from 
the uncertainty in the measurement of \q/p\ and 0. The constant g^^-^ can only be determined empirically. As 
for g^-^, the contribution of the El channel is small for both the decay rate and the asymmetry, and hence the 
uncertainty stemming from it is not significant. 

Despite these large uncertainties, it is possible to set an upper bound on the size of the asymmetry one can 
expect in these decay channels. In the analysis to follow we shall show how the upper limit on the asymmetry 
varies with the uncertainty in the magnitude of each parameter. As for the phases, since there is no concrete way 
of constraining their range, to get the upper limit on the asymmetry, we shall assume maximal phase contribution 
to the asymmetry which in mathematical terms means, given 

= Ai cos Oi + A2 cos 62 (62) 

TT 

©1 = arg(.gMi5BR) = i'^± + So - Si - -) mod tt (63) 
02 = arg(.gMi5Ei) = (*± - ^) mod tt (64) 

we shall assume Qi ~ O/tt and O2 ~ O/tt. It should also be noted that we are not varying the parameter gp due 
to our lack of theoretical control over it. Suffice it to say that a smaller than estimated magnitude of gp will raise 
the maximum value of the asymmetry from 0(1%) to ~ 10%, the phase contribution remaining maximal. 

In Fig. 3 we plot the variation of the asymmetry over the range of possible values of g^^i and g^p. For electron 
pairs in the final state the asymmetry can be as large as 2%. However, it should be noted relatively large asymmetry 
is found for relatively smaller values of g^p^ and large values of .9^1 ■ This is contrary to the plot in Fig. 3 that refers 
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9 Dl ^ h+h-l+l- 




\9m\ 



Figure 3: Contour plot of the CP asymmetry parameter A!^ with pious in the final state in the space of the 
couplings parametrizing the Ml and bremsstrahlung contributions. Results for both electrons and muons in the 
final state are shown above. The dashed line marks the upper limit on gsR when assuming there is no direct CP 
violation. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



\9m I I^Mi I 



Figure 4: Contour plot of the CP asymmetry parameter Ai^ with kaons in the final state in the space of the 
couplings parametrizing the Ml and bremsstrahlung contributions. Results for both electrons and muons in the 
final state are shown above. The dashed line marks the upper limit on (7br when assuming there is no direct CP 
violation. 
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to muon pairs in the final state where the larger asymmetry occurs in the region of larger ggp and 5mi- Also, 
the maximum asymmetry in the decay channel with muons in the final state is smaller by orders of magnitude as 
compared to the ones with electrons in the final state. These, in general, seem to hold good for kaon pairs in the 
final state as seen in Fig. 4. 

Both the fact that the asymmetry is much lower in the case of muonic final states and its dependence on g^^^ and 
5br '^^^ understood from the finiteness of the muon mass as compared to the electron mass. The bremsstrahlung 
component of the amplitude is sensitive to the mass of the final state leptons and there is a natural cutoff for this 
contribution from the muonic mass which serves to naturally curtail the infrared sensitivity of this amplitude at 
low leptonic invariant mass. However this is absent for the electrons and large infrared contributions are seen in 
the bremsstrahlung component and its interference with other amplitudes. 

The numerator in the asymmetry depends linearly on the bremsstrahlung contribution, but the denominator 
depends quadratically on it. Hence for the electrons the asymmetry peaks and lower values of before the 
quadratic dependence in the denominator takes over and reduces the asymmetry. On the other hand, for muonic 
final states this does not happen as the bremsstrahlung contribution faces a natural infrared cutoff at the dimuon 
mass. 

From Fig. 3 and Fig. 4 it is evident that assuming that there is no source of direct CP violation in charm dynamics 
does not affect the possible size of the asymmetry to a large extent. While for the case of electrons in the final 
state, the region of maximum asymmetry lies well within the upper bound set by the dashed line, for muonic final 
states regions of larger asymmetries are not completely excluded. 

9.6 e+e- [h+h-][l+l-h+h-], [Ks4>][l+l-h+h-] 

While even the existence of e'*"e~ — >■ [D'^ D'^]c=- [h'^h~][l'^l~h^h~] is a signal for CP violation the same is 
not true for e+e" — )• [D°l5°]c=_ — )■ [K s(t)][l'^ h'^ h^]- However, a detailed study of correlations in both the final 
states is necessary. Since the final states are common to 13° and Z)", the integrated partial width of the correlated 
states is given by Eq.38 

10 Other four body decays with a lepton pair in the final state 

Within the SM and the KM ansatz, CP violation can manifest itself in Cabibbo suppressed decays of charmed 
mesons. However, as mentioned earlier, there are some subtle exceptions. CP violation can exist in the DCS decay 
_D" /v+TT^ vs. L*" — > A'^7r+ and the apparently Cabibbo allowed decays Kstt^- 

D° K+7T-1+1- vs. D° K-7T+1+1- 

CP violation can manifest itself in Z?" — >■ K~^Tr~ vs. — )• K'n'^ through the interference of AC = 1 
dynamics with time dependent CP asymmetry as mention earlier. An enhanced asymmetry can possibly be 
observed as in the case of U"™* — > h^h^l^l^ , as discussed above, through the measurement of a T odd 
correlation although the same price has to be paid: a much reduced branching fraction. 

D+ Ksh+l+l- vs. D- Ksh-l+l- 

For KsTT^, an apparently Cabibbo allowed mode, the final state can be accessed through a DCS 

reaction through K^—K^ oscillations and the two amplitudes interfere, opening up possibilities of CP violation 
even within the SM and, more importantly, allows for possible ND contribution as mentioned earlier. This 
can also lead to a T odd correlation in the comparison of — > KsTr^l^l~ ■ 

In the case of — > KsK^, a Cabibbo suppressed channel CP violation is manifested only through the final 
state and is quantified by 3fi(e_fs-). This effect will be reflected in the comparison of KsK^l^l^ . ND 

cannot leave any impact here. 

11 D"^"* hli^ and D hvi? 
11.1 l-ph+ vs. D° l+uh- 

The rate of both l^vh^ and 13° l^vh~ are given by the same factor e"'"*, but suppressed by the tiny 

4 4 

factor (a;|, + y'j^)/2. Their asymmetry is given by ^4~^4 independent of their time of decays. Yet the experimental 
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bound on \q/p\ allows the range from about 0.73 to 1.1 and allows the to range from -0.56 to 0.19 - i.e., 
a possible large asymmetry. Furthermore non- ad-hoc models of ND - like LHT - allow a significant part of this 
range, in particular if it produces a sizable, but not leading part of the observed value of a:;^,. The same is true for 
models with non-minimal Higgs sectors. 

As stated before this asymmetry does not depend on the time of decay of the neutral charm meson. Yet two 
comments should be mentioned here: 

• The data, on which the observation of such an asymmetry are based, should show the 'known' time dependance 
for both and D'^ transitions. 

• If ND produces 'wrong sign' decay through AC = 1 dynamics, it would generate different time dependance 
of D° l-X vs. 5° ^ l+X. 

It is crucial to determine whether the neutral meson was a or . The usual flavour tagging is done by 
D*+ D"7r+ vs. D*~ D^TT-. Another way would be to employ B- D^oi*^^ and B+ D^oi*^^ to tag 
vs. l)" as discussed before. 

11.2 Asymmetries in Dl,s l^uh^ 

One can consider a search for a CP asymmetry in semi-leptonic decays of D mesons in analogue with kaon decays, 
where one has measured the predicted difference in — > l^vn^ . However the analogy is at best only qualitative, 
not quantitative. For neutral kaons one has large oscillations due to AMk/^k and a huge difference in life times, 
while Hqk/pkI — 1| < 0.005; for neutral D mesons we have small values of x^,, y^^ in the range 0.005 - 0.01, 
while Iqd/pdI could differ from unity by 10% or more. In e+e~ annihilation close the charm threshold one can 
study e+e~ — ^ DgD^ denoting the two mass eigenstates with short and long lifetimes and make use of quantum 
correlations ~ as one had for e+e^ — ?> 0(1020) — > KsK^. Of course one has to consider the hurdle of small 
lifetime differences in the neutral D system but can reap benefits from the possible large CP impurity in the mass 
eigenstates. 

11.3 D hui? 

In a previous work [19] we have discussed rare charm decays like D — > l^l~ X and found that LD dynamics have 
a strong tendency to overwhelm signals from ND, unless one finds CP asymmetries there. Yet the situation is 
different for D — > XvD, since LD dynamics cannot contribute there in an appreciable way. SM dynamics can 
generate D — > XvD, yet with branching ratios well beyond possible experimental reach. A simplified calculation 
using the Inami-Lim functions [20] leads to the numbers [21]: 

BRsd(£»+ ^ A:„i/P) ~ 1.2 X 10~l^ BRsDp" ^ ~ 5 X 10-1'^ (65) 

However, to treat this decay channel properly one has to take into account the fact that the quarks running 
in the loops of the penguins and the boxes are essentially light quarks and cannot be treated like the top quark 
loop as is done in the analogous B or K decay channel. To reduce theoretical uncertainties one would have to do 
a complete two-loop renormalization group analysis [22] as is done for the light quark loops in the analogous K 
decay channel which, considering the inaccessibility of these channel to experiments, is not a quite fruitful exercise 
at present. 

On the other hand, since these decay modes are highly suppressed in the SM and get negligible contribution 
from LD dynamics within the SM, they are ripe for ND intervention. Our analysis of these channels within a 
LHT-like^ ND reveals that enhancements of 0(10'^) — 0(10"*) is possible in a large part of the considered parameter 
space'^ as can be seen in Fig. 5. Even with such enhancements, these decay modes will remain beyond realistic 
experimental reach for quite some time. 

Within the SM, there can be no CP violation in D Xvv. LHT-likc models can do very little to change this 
as direct CP violation remains extremely tiny and mostly nonexistent in such ND. 

^For the details of the ND c.f. sect. 13 

description of the parameter space can be found in [19,35] 
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12 Comments on £) — >■ hih2hsh4 



Constrained by B and K Physics Not Constrained by B and K Pliysics 




Logio[BR(D ^ X„vv)] Logio[BR(D ^ X,vv)] 

(a) ruHL = iOOGeV 



Constrained by B and K Physics Not Constrained by B and K Physics 

t 




Logio[BR(D ^ X,vv)] Logio[BR(D ^ X,vv)] 

(b) ruHL = llOOGeV 

Figure 5: Histograms of BR(I? X^j^vv) including LHT contributions for mirror lepton mass of (a) niHL — 40006^ 
and (b) ruHL = llOOGey 

Furthermore the mere existence of Dl tt^vv would already show CP violation. At the leading order this 
process can occur only through penguins and box contributions, both of which would be CP violating 
channels. The argument is clearer when one looks at the effective hamiltonian 

Under a CP transformation (j)T^o — >• — 0^o, d^(j)jj^ — >■ — 9^</>_d^ and 1/7^1' — > —v^^^v making "He// CP odd. Since 
there are no higher order contributions to this channel (unlike — > Tr^e+e^), in the absence of CP violation this 
channel is forbidden. Hence, any finite measurement of this decay mode is a clear signal of CP violation. 

12 Comments on D —Y hxh2h^h4^ 

Obviously there are theoretical and experimental challenges in treating four-body final states. However the analyses 
of four-body decays offer advantages too, namely there are several obscrvablcs of CP asymmetries that can (a) offer 
checks on measurements and (b) can give complementary information about the underlining dynamics. We give 
some examples in non-leptonic D decays. 

Let us consider "tt+tt". As usual, a time dependent CP asymmetry can be studied here by 

comparing D*^ — )■ K^K^ t:^tt~ vs. D*^ — )■ K^K^tt^tt^ . However there are other probes that can be employed 



20 



13 A non-ad-hoc scenario of New Dynamics - the Little Higgs Models with T Parity 



very effectively here. The T odd correlation is defined as: 



Ct = {{pk+ X X p^-)) 



(67) 



which under time reversal transforms as: Ct — > —Ct- However, Ct 7^ does not conclusively establish T violation 
as it can be driven by FSI through T invariant dynamics, time reversal being an anti-unitary operator. The effects 
of this FSI intervention can be resolved by measuring the conjugate T odd correlation: 



as Ct 7^ —Ct unambiguously proves the existence of CP violation. 

The effects of ND can appear in this asymmetry without having left a large impact in the partial decay rates 
themselves. Arguments placed in sect. 6 also govern the size of the asymmetry here. A comprehensive study of this 
asymmetry with LHT-likc ND is not presented in this paper and is left as a future task. Some interesting comments 
on such four-body hadronic final states can also be found in [23]. 



13 A non-ad-hoc scenario of New Dynamics — the Little Higgs Models 
with T Parity 



Flavour physics has often served as testing grounds for any ND and has played a majors role in the christening 
of different models as plausible candidates for dedicated searches in some confines of their allowed parameter 
spaces. The Little Higgs Models [24-27] are no different. These models were originally formulated to address the 
hierarchy problem through extensions of the SM global and gauge groups and the judicious arrangement of their 
couplings along with collective symmetry breaking, something which is complimentary to how Supersymmetric 
(SUSY) models address the same problem, through transformations of spin statistics. However, like in SUSY, a Z2 
symmetry, called T parity [28,29], had to be postulated to limit contributions that could break the SM custodial 
SU{2) symmetry which, in turn, would allow the scale of ND to be brought down to around 1 TeV . 

One of these class of models found favour amongst theorists for its minimal structure: the Littlest Higgs 
Model [30] with T parity, having a relatively small global and gauge structure and only a few free parameters. 
One implementation of T parity in this model necessitates the introduction of a set of "mirror" fcrmions [31-33]. 
Although the model was built to address the hierarchy problem, to make its implementation consistent required 
the model to be naturally non-minimal flavour violating, a feature which is much sought after by flavour physicists. 
This feature then opens up avenues for possibly interesting ND features in rare decays and CP asymmetries in K, 
B [34] and D [4, 13, 19,35] decays, a large part of which has been extensively studied and exists in the literature. 

What makes the flavour structure of this model particularly well studied is that a significant space of the same 
favourable parameter sets have been studied m. K , B and D systems generating significant correlations between 
such ND contributions to oscillations, rare decays and CP violating parameters in all three mcsonic systems. These 
studies not only bring together the "up" and the "down" type in the SM, but also the "Up" and the "Down" type 
in the mirror fermion sector. 

For the D meson system, it has been shown that LHT can account for all the strength of the oscillations in 
the neutral system [13]. It has also been shown that LHT can generate very little direct CP violation in the 
inclusive leptonic channel and two body final state hadronic channel [4,19,35] with SM long distance contributions 
overshadowing contributions from LHT to decay rates in the rare decay modes. However, the story is very different 
for indirect CP violation where LHT can leave its mark through its contribution to the oscillation and mixing 
parameters. The philosophy which governs our view of ND intervention in the two, three and four body decay 
modes that we discuss in this work is based on the fact that in the light of the uncertainty that surrounds the SM 
contribution to indirect CP violation in the D meson system, LHT-like [4, 19,35] ND can step in and make some 
or all of the contribution possible. 

What has been important in all our previous analysis on such ND in D systems, and is also pertinent to the 
current analysis, is the fact that we do not see LHT in itself as being the solution that nature chose, but rather, 
when stripped of the detailed constructs of the model there lies in it a very interesting flavour structure that might 
be generalized by certain principles elucidated in [4, 19,35]. Such flavour structure can be generalized to a few other 
models that have already been postulated or can appear in other flavour structures that might be formulated in 
the future. As our analysis is based on general constraints and extensions of SM like flavour structures, we feel 



Ct = {{pk- X (pV-) X p^+)) 



(68) 
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15 Summary 



that our conclusions go beyond LHT itself and encompasses any model that might share with it similar flavour 
dynamics. 

14 Comments on ND with non-minimal Higgs Sectors 

Having considered a class of models that can leave large impact on indirect CP violation while leaving almost no 
footprint on direct CP violation, it is natural to consider another class of models that can affect both. One such 
class of models can be built by introducing non- minimal Higgs sectors in the ND. Such models typically have extra 
scalar and pseudoscalar degrees of freedom and CP violation can be introduced through complex Yukawa couplings 
and/or mixing between the scalars and pseudoscalar higgs. 

Depending on the details of the Higgs sector, i.e., the number of Higgs doublets and the nature of the Yukawa 
couplings, different implementations of a non-minimal Higgs sector have different contributions to both FCNCs 
and CP violations in the quark sector. Varying prescriptions can be used to reduce the contribution of these new 
Higgs sectors to FCNCs to make them compatible with electroweak precision measurements and flavour physics 
constraints. Most common prescriptions lead to a MFV structure. However, it is a common misconception that in 
such MFV constructs there can be no new sources of CP violation'*. While dangerous contributions to FCNCs can 
be avoided by the prohibition of new sources of tree level contributions to FCNCs, CP violation beyond the CKM 
construction can arise from new flavour-blind phases. It is also possible to arrange for some simple symmetries in 
some models with multiple Higgs doublets such that FCNC contributions to the down type quarks is limited to 
within the experimental bounds while allowing significant contributions to the up type quarks [37]. 

Unlike, the LHT-like models considered above, models with non-minimal Higgs sectors can contribute to both 
direct and indirect CP violations. While the presence of a charged Higgs along with possible complex phases in 
the Yukawa coupling can induce tree level direct CP violation, the running of these scalar/pseudoscalar degrees of 
freedoms in loops can induce indirect CP violation. For the case of the neutral D meson such sources of indirect CP 
violation can saturate the experimental bounds on the oscillation parameters. Moreover, charged Higgs currents 
can leave tree level contribution to direct CP violation in tt+tt^ and more so in the DCS decay — ?> K^ti^ 

where SM contribution is nonexistent. As we have seen above, CP violation in — >■ h~^h~l~^l~ is driven by the 
CP violation is Z?" — ^ h~^h~ and hence enhancement in the latter can easily manifest itself as an enhancement in 
the former. A detailed study of such effects in charm physics is left as a future task. 

15 Summary 

Experimental inaccessibility of almost all rare decay channels and asymmetry parameters have in the past led to 
a lot of skepticism about the viability of both theoretical and experimental study of charm in the same glory as 
granted to beauty and strange physics. However, that indirect CP violation in charm is ripe for ND intervention 
has been increasing finding favour amongst flavour physicists. This has been fueled primarily by the conclusive 
observation of charm oscillations along with a realization that SM in all probabilities cannot account for the same. 

Armed with this knowledge, we have proceeded to establish a generalized flavour structure, the LHT-like flavour 
structure, well tested in B, K and D physics, as a contender for showing effects in two-, three- and four-body decay 
channels in charm decays. Although, there is little hope that such ND can enhance direct CP violation, indirect 
CP violation can show up with very unique features, as a very small effect but with sizable CP impurities in the 
neutral D mass eigenstates. We have also proceeded to comment on a different class of models with non- minimal 
Higgs sectors which, unlike the LHT-like models, can bring in both direct and indirect CP violation. What effects 
they can leave on the kind of charm physics discussed in this work is yet to be established concretely. However, it 
goes without saying that both kinds of ND discussed above can saturate the current experimental bounds set on 
CP violation in charm. 

As we have emphasized in our previous works on rare charm decays within a LHT-like flavour structure, our 
conclusions arc more general than the model that we have tested it in. The same is true for the current work. 
In addition, decays with higher multiplicity in their final state open up access to more experimental observables 
although they are theoretically more complicated to deal with than decays with lower multiplicity in the final state. 
Even with the advent of Super B factories and the LHC where copious numbers of D mesons are produced, tiny 

*For an intriguing counter example see [36] 
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branching fractions and large backgrounds, specially for the case of hadron colliders remain obstacles that will have 
to be tackled in the years to come. 

Despite all the effort that is necessary to glean experimentally observable effects in the D meson system, it is 
already showing promises of unique effects that have hitherto been unseen in both the B and K systems which 
have been studied extremely well both theoretically and experimentally. Along with this, the fact that charm lies 
ripe for the intervention of ND should be sufficient to keep some intrigued souls hard at its heels. 
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